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I. INTRODUCTION 


The decomposition of the Laplacian operator, 


^ = F 3 I? ( r2 sF ) + F 3 (jet? t^ sln 9 h + ' 


2_ 



A 


where the coefficient of the second term is proportional to the 
square of the angular momentum operator, is the basic relation 
between energy and angular momentum in the quantum mechanics of 
the one body problem (or the relative motion of two particles). 

When acting on a wave function which is an eigenfunction of total 
angular momentum t, the Laplacian simplifies to 

V 2 = f- fr 2 |-V i- l) 
r 2 ar V ar J r 2 

in which form it is clear that the effect of this decomposition is "l 
to reduce the SchriJdinger equation from a three dimensional partial 1 
differential equation to a one dimensional (ordinary) differential 
equation. As such this relation is of fundamental mathematical 
importance and is familiar to everyone who utilizes quantum 
mechanics at all. 

The analogous procedure when more than one particle is 


GPO PRICE $_ 

p c CTI DDIPCIC\ 

j 

N65-35271 


Lorll rrilUfcvo) 

1 

3 cm 

! (ACCESSION NUMBER) “ 

i 

(THRU) 

/ 

Hard copy (HC) 

(PAGES) 

i r7M\ - £T/5Z)5 

(CODE) 

Mir.rnfir.hp (MF) 

3b 

(NASA CR OR TM* OR AD NUMBER) 

(CATEGORY) 


1 • 



- 2 - 


involved, in particular two identical particles in an external 
force field, although known, is not as well known, nor is it as 
well developed.! When the external field is that of a fixed nucleus, 
the wave function is expanded in eigenfunctions of the total angular 
momentum of the two particles multiplied by functions of the three 
remaining independent variables. The total angular momentum eigen- 
functions are functions of the three Euler angles only. These angles 
are not unique, but in some way they must describe the orientation 
of the instantaneous plane formed by the two particles and the center 
of coordinates (nucleus) in space. The remaining three coordinates 
then describe the positions of the particles in this plane, and the 
functions of these variables are the generalized radial functions. 
Hylleraas ' original papers'*" in effect contained the reduced or 
radial equations for total S states in terms of the residual coor- 
dinates ri, r 2 ,ri 2 « In this case, the total orbital angular 
momentum is a constant function, and hence the reduction of a six- 
/ dimensional to a three-dimensional partial differential equation is 
independent of how one defines the Euler angles. 

2 3 

The standard treatment of the general problem is due to Breit, 1 
He used the Euler angles that Hylleraas 1 originally introduced: 
namely the two spherical angles of one of the particles in the space 
fixed coordinate system and a second azimuthal angle between the 
ri - z plane and the r 2 -ri plane. Breit' s remaining coordinates 
were chosen as ri, r 2 , and 0 i2 , the latter being the angle between 



r x and r 2 . If, however, one wants to describe two electron atoms or 
ions in the approximation that the nuclei are fixed, then one has 
an additional requirement of which there is no analogue in the one- 
body problem. And that is the Pauli principal: the requirement 

that the spatial function be either symmetric or antisymmetric under 
the exchange of the particle coordinates. It is clear that the 
Hylleraas-Breit choice of Euler angles (which we hereinafter refer 
to as the Hylleraas-Breit angles), being quite unsymmetrical with 

if 

respect to the two particles, is not optimum in this respect . In 
fact the construction of linear comb ina tions of an gular momentum 
functions with the appropriate exchange properties is a very difficult 
task which depends not only on the Euler angles but on 0 12 as well. 

p 

It is not surprising, therefore, that Breit's original work was 

limited to P-states, and work thereafter has always been limited to 

5 

specific angular momentum states . 

However, a treatment by Holmberg^ using a symmetrical choice 
of Euler angles (which we shall call Hdmberg's angles) has in the 
interim been carried out. With these angles the description of 
exchange as well as parity (which latter property is also simply 
describably with the H-B angles) is simple (although these proper- 
ties are only alluded to in Holmberg's paper'’ ) . One of j the pur- 
poses of the present paper is to examine these properties and relate 
them more clearly to the construction of the total wave function 
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an<i further clarify other aspects of Holmberg's important paper. 

Principally, however, we shall: derive the general radial equations 

for arbitrary angular momentum (JL ) for the case of two identical 

particles in an external field. Holmberg's treatment applies to 

three particles of the same mass. 

It is clear that in treating three particles of equal mass 

Holmberg had in mind the three nucleon problem whereas we are 

interested in two-electron atoms, ions, and diatomic molecules. 

The application of this formalism to two electron atoms and ions 

is clear, and the decomposition amounts to a rigorous reduction of 

.the SchrdS dinger equation. It should only be remarked here that the 

scattering of electrons from one electron atoms and ions is also a 

special class of these problems. J We have therefore worked out the 

connection between the boundary conditions for electron-atom 

scattering and Holmberg's angles (Section VIII). j 

Inasmuch as Holmberg's paper refers to the three nucleon 

problem, reference should also be made to the papers of Derrick and 

7 

Blatt . These papers deal much more realistically with the three 
nucleon problem in that full account is taken of an internucleon 
potential which is considerably more complicated than a central 
potential. As regards the actual choice of coordinates. Derrick and 
Blatt define axes along the moments of inertia of the three -body 
system. As such they will depend on the lengths of the interparticle 
distances and therefore are quite different from Holmberg's angles 
(Section II). 
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II. HOLMBERG'S ANGLES 


Figure 1 contains a perspective drawing of Holniberg's angles 
which define the particle plane with respect to the space fixed x, 
y, and z axes. The rotated axes x', y', z' are then defined by 


2 ' 


X' 


l£i* $z\ 


z x ■£' 

MV -V* 


z -X 


2' 


( 1 ) 

( 2 ) 


= £' X x' 


( 3 ) 


/ 


The Holniberg Euler angles are then 

0 = angle between z and z.' 

<§> = angle between x' and x 
yy = angle between x_' and (r 2 -ri) 

The ranges and planes of these angles are: 

0 ^ 0 S- ^ in z-z' plane 

0 4 $ 2 "h" in x-y plane 

0 5 H $ 2 "W inx'-y' plane 

As is usual a cap on a vector is used to represent a unit vector 

A A A 

in the given direction. In particular i, j, k are the three 


W 

( 5 ) 

( 6 ) 


unit vectors along the (space fixed) x, y } and z axes respectively. 



A ^ ^ A ^ 

and thus are synonymous with x, y , and z_. Similarly V , j ', 

A A A A 

k' and x', y , zj are identical. 

It is clear from the figure that x, "being in the x-y plane 
has components v 

&' = C # 4- 3 S ~ ^ 

Since x 1 is perpendicular to the z-z ' plane, it is perpendicular 
to every line in that plane going through the origin. This in- 
cludes specifically the line of intersection of the z-z' plane with 
the x-y plane. However the azimuthal angle of that intersecting 
line is the azimuth of itself, and since x ' has azimuth <£ , z ' 
hag azimuth ~~ + ^ (cf. figure l). The polar angle of z' is 
clearly ©, therefore we have the important relation: 

A ^ 

£ _ l Si. e s~ $ - i s~ e • Cc, <£ + k g.s e 

The relations between Holmberg's angles and the spherical angles 
of the individual particles are obtained by substituting Eqs. ( 7 ) 
and (8) into the Ihs of Eqs. (l) and (2) and using the ordina r y 
decomposition of rj. and r^ in the rhs : 

A A ^ A 

/i 2 ^ 5^ ^oi ^ Sii T? Sv«t 

r r - t + \ Sw ^ s- • 


One obtains 
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Sv» 6,^ Coi8 = S~ S _ «p ( ) 


(9) 


S“*®ix i>i& Si$ - Sw 1 ^, SiK<fJ <-o s tf x - C.* 2f, Si lfc* Si. <p^ 


(10) 


^*-^ 11 . S«i 6 Co5$ - Si i?, Co^ — G.t j) ( S.'« i? a CoiCp 

2 Q IX ) CaW = Si 2s a e. »(<?*-$ ) - Cdo^-^) 

* SH^O SwV ^ (Cc^ -c. s ^ ( ) Si0 

-s^* ( :*/-(<£- £ ) + Si^ Si (c^-4) 

(^> e i-v = Coi Cos Z\ -*- SiT>, S~- 


( 11 ) 

( 12 ) 

(13) 

(14) 


The latter relation is, of course, the well known expansion for 
the angle between two vectors. 

It is also of interest to give the vectors r x and r 2 in the 
particle plane: (primed co-ordinate system): 

r = HSiiw- ie.,1 - \ c., ly- i«n) 

A ^ 

S - i + ie a ) - £ Cc S (V + JLfc fl ) 

A r, A/ 

Z- ~ 0. Sw © +• t CoS 0 


(15) 

( 16 ) 


( 17 ) 
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^ = C_ $ +■ Cos, 6 < J' — k 0 Cos $ 

- C Cos ^ — J '6 S'- ^ ^ K S«« 0 S<i« $ 


(18) 

( 19 ) 


The following relations which are also very useful can now simply 
he derived by computing (£i • z) and (r 2 * z) in the primed system. 


Cot dj ■» — Sa» 0 C.J ('Jf — ®l«.) 


(20a) 


Cot a — Sw 6 C -1 I'j' + (20b) 



S%*» 

II 

Sr 

•a 

o 

Cos ^ - C.jS Cot ^^<1-) 

(21a) 


Si m ^ 

C*J ^ , 

Cm$sM¥+ - c.»9 Sv-<| c«o»(u + ^6,,.) 

(21b) 

- 


S~ <f ( , 

S*~ l u - C».s 6 Cos^ Cos(sjj- 

(22a) 

- 


- 

s.~ (.m: -v - c< ‘ s ® c * s< % Co ' ^ 

(22b) 


III. PROPERTIES UNDER PARITY AND EXCHANGE 


The operation of parity corresponds to the simultaneous in- 
version of both particles; coordinates: r x -^ -ri, r^ ->-r 2 . It 
can be seen from figure 1 that this places ri and r 2 facing the 
opposite direction, but the cross product and hence z' will not 
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change as a result of this operation. Thus the z-z’ plane will 
not change and x' will not change. On the other hand (r P -ri) 
goes into the negative of itself, so that V gets increased by ir . 

In other words under parity 

0-^0 

9 -+ £ (23) 

<y/ + IT 

Exchange corresponds to the transformation ri% £2* From the 
analytical definitions £' and Eqs. (l) and (2), the new primed axes 
will go into the negative of themselves. Also (r 2 -rj.) goes into 
negative of itself. Clearly the inversion of the z' axis corres- 
ponds to the transformation 0 tt - 9. Noting that $ is the 

angle in the x-y plane and measured as positive with respect 
to the z axis, which is fixed, we see that f § . The 

simultaneous inversion of x' and (r 2 -ri) means that the modulus 
of the angle y remains the same. However, since "f is an angle 
in the x'-y' plane, it is measured as positive with respect to the 
z ' axis. Since the latter goes into the negative, itself , it 
becomes clear that \ -+ xt - y • Thus we have under exchange 

9 -» TT- 9 
^ (^> -v "N" 


(2k) 
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The significance of these transformations relates to the 
transformation properties of the vector spherical harmonics under 
the same operation. These functions, which are the eigenfunctions 
of the angular momentum (next section), are the basic functions in 
terms of which the complete wave function is expanded. They can 
be written 

. .(■«#«•*<¥) -k 

* ja-ut+l? e, JL.ie) 

4-TT 


(25) 


where the normalization has been so chosen that the function is 

y*,k 

identical with what is given in section IV and the c(^9) agree with 

8 

those given by Wigner . Only the dependence on 9 is non- trivial; 


for m j/ k J/ 0 


at**- 7 " 


m - 

SC. 1- S 

x 


j m 'rV> - 1 a- c.» i 3 

- 4 (TTk)« 

* - k ♦ » , — Van i. b) 


(26) 


F(a,b; c; z) is the hypergeometric function in the notation of 

o 

Magnus and Oberhettinger . The important property of d t ^©) , proved 

Q 

in Wigner 's book , is: 


d f C*-e) = t-\) 


(27) 


Letting P and £ represent parity and exchanges we have 
from (23) and (24): 

P 2>t (9, $ > V) = 

v)= ^(ir-e^ + TT, 2TT-V) 
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which using (27) reduce to 

K 


y* , vl 


pa, ie, <f , Y) C-' )* 3l t U ,<t .V) 

£u V) 


C 

* C-0 % Co ,$>V) 


The simplicity of Eq. (29) is the essential feature which 
recommends Holmberg ' s angles to the description of the two electron 
problem. 

IV. ANGULAR MOMENTUM 

The components of the total angular momentum are readily 
expressed in terms of the particles' spherical angles. Thus, for 
example 


-k i ♦ Cot-0, C «-><*> > + cp 2 

*- 1 — ai 




2 

fl9 t 


+• Cofc l? x CaSCfi^ 

dC fu 


The particles' angles "6i, & 2 > *¥ 2 vla (9) ~ (14) 

are implicit functions of the four angles ©, $ » V > ®i 2 * 

Thus the problem of finding in these angles is a straightforward 
problem of partial differentiation. We can write 


(28) 

(29) 


(30) 


/(- * S^-d> 4- 0o4 &, Cgj <? t CoH^C®a^^2f' 




?<e, 


*1*. (31) 


where 
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and f' can be anyone of the angles 9, $ , V , or 9ia* Using 

then Eqs. (9) - (14), one finds that the following relations fall 
out quite easily: 

A e r O 

A 0 — — Cos <|> 


(32) 

(33) 



Svw ^ CoS 0 


(34) 


A 


h: « 


<£ 

s^e 


(35) 


Thus 

M x * [ Cos $ - S«<$Coi& e> + 

L *e -3$ su,e 

One of course, proceed in a completely analogous way to 
get the remaining components of the angular momentum, however, let 
us note from Eqs. (20) that 

<£#1 =. * o 

and from (10) and (ll) 

yf, _ ils 8 | 


l 1 
w J 


(36 


(38) 



Since 


a 


I 




( Vti . i 


+ . 5 ^ 

*£ at^ y 


substitution of (37) and (38) 


yields 


1 . 2l + ^ 

However, since 

K f iL + 2 . ) 


we therefore have the z -component of M: 

M- . * 

* c 

The remaining component of the angular momentum may be derived 
from the comnutatlon relation H 

Straightforward substitution yields: 

— fe T" Su. $ a. -i- Cos <£ c«t © — ops j t *a_ ~J 

<r t L ^ ® Sv« © ^y* •* 

These relations sure independent of ©i 2 corresponding to the 

statement that the angular momentum only depends on the (three) 

Euler angles 9, . The forms of the three operators is the 

3 

same as one gets with the Hylleraas-Breit angles . The square of 
the angular momentum is likewise the same. One finds directly 
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from the sum of the squares that 


M - — t) l" 6 ■+■ Cn fc 0 ) £. i. +■ J_ 

L -df x '6$'- s^0 yy*- 

— * ^°t_P £ 0 £. 9 <) ~j 

^uv, e ■zi'b't J 

The vector spherical harmonics, which have been given for 
the restricted range m ^ k 7/ 0 in Eq. (25 )> are the similtaneous 
eigenfunctions of with eigenvalue % x {,(. l -h) and ^ ™ 

.***. k 


M (e, $ , V) = V IU*0 3 c ('e , § , 'It') 

My 3 £ , $ 90 - t, >vi 4 (8 , <§ , ty) 


( 41 ) 


(42) 

(45) 


They are given in a completely general, normalized form in Pauling 


and Wilson' 
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|k-H / k + W L(vn$+k<^) 

Coi > e e 


, FC-e + i/a-i, * * ) 


(44) 


where 


f& = i Jk + m + k - m + 2 


(45) 


and the normalization constant is 


N 


/hK 


+-L lk + »"l + ^ I K -»»))! ({. - l\k + -*»-\ -v L | ic - >«0 1 

- ^|*->«|) • "* -t|K- »*»()! i) 


1_ 

-N 2 - 


(46) 
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In addition to the usual magnetic quantum nurrber m, the 
vector spherical harmonics depend on the quantum number k, an in- 
teger whose range of values is the same as m: - t m, k t 

The physical significance of k derives from the fact that the <3^ 
are the eigenfunctions of the spherical top (for which (42) is the 
Schrfldinger equation), and k is the angular momentum quantum number 
about the body-fixed axis of rotation. With regard to the 
applications that we contemplate here, k can be considered a degen- 
eracy label which must be adjusted such that other requirements are 
fulfilled. 

V. CONSTRUCTION OF THE TOTAL WAVE FUNCTION 

We shall confine ourselves here strictly to the atomic problem vhich 
implies that the potential energy as well as the kinetic energy 
commutes with the total angular momentum. In this case the total 
wave function for a given Jt must be a linear superposition of the 

Cl’”./* 

degenerate ^ . In addition, m will be fixed for a given mag- 

netic substate and the "radial" equations will be independent of m 
(cf. Appendix II). 

Considering, for the moment, the residual coordinates as ri, 
r 2 , ©i 2 » we can therefore expand the total wave function in the 
form: 
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The parity operations, Eq. (25), only effects the Euler angles, 

k 

and from Eq. (28) it only multiplies the ^ function by (-1) . 
Therefore by restricting the sum to even and odd values of k, we 
guarantee that the superpositions have even and odd parity respec- 
tively: 

c 5 . *•> = 1 _ 1/ v n .e,o 31,(6. <J ■ 

k 


Octet 




'1) 


1 J,V. .e ,e,») 

K elV 


(48a) 

(48b) 


where the double prime on the summation emphasizes that the 6um 
goes over every other value of k. 

In deriving the radial equations (next section) we shall ex- 
ploit the invariance of the radial equations with respect to m, 
by choosing m = 0. When the Hamiltonian is written in terms of 
the Euler angles and the remaining variables, there will occur Serins 


involving 1. and iL . By virtue of m = 0, the former terms 
!>£ W 

vanish, but the latter terms would bring down the imaginary co- 
efficient ik. In order to avoid complex equations, it is therefore 


convenient to construct real angular momentum functions. Let 




- 


x (3 1 ; 

J5 




2T*) £7. 


o. % 


+1 Cot Xt|/ dL. (S) 


xir 




21T 


• to) 


(49a) 


(49b) 
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for x y 0, where 


r = im 


(50) 


and for Y = 0 define 


2 °;* *7 ■ 


(49c) 


This then constitutes a set of real, orthonormal vector spherical 
harmonics (note = d^ ^). These real vector spherical har- 

r 

monies are still eigenfunctions of parity with eigenvalue (-1) 

The property of exchange is a mite more complicated than parity 
in the sense that it affects not only the Euler angles, but the 
residual coordinates as well. The beauty and importance of Holm- 
ber's angles, however, is that there is no mixing, and independent of 
whether we consider the residual variables ri, r 2 , 612 or ri, r 2 , ri 2 , 
the effect of exchange on the residual coordinates is simply r x -*r 2 . 
Finally then, if we construct y+ 

♦ r ^ * 

\ j/ fo cs .S) - Z y) 


+ ,e.o 3/®,# , ^ >] > 


(51) 


the operation of exchange on this sum then gives, with the use of (29), • 
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*i io (£•.?.) = ( t , 5 ) 

r X+ ( X+ 

= I [4/t, ,£,&,) C-0 (o.i.v) 

*- *' 

+ h • £ - e ”) ') *<V e ' ^ 

Thus if 

f x+ 

,0,0 = ± C-o • f x (Y,, Vi.eu) 


f+4 tf- 

“7 . 6 ta ) = ± C-i) ^ ( l^i • / ©i, ) f 


(52a) 

(52b) 


the function \i/ / of Eq. (51) is a real, space symmetric (upper sign) 
to 

or space antisymmetric (lower sign), eigenfunction of M 2 and M cor- 
responding to the quantum numbers X an( l ® with m => 0. The space 
symmetric and antisymmetric solutions correspond to singlet and trip- 
let spin states respectively. Furthermore the restriction to m = 0 
is sufficient for deriving the radial equations. 

We have shown that the m = 0 function can be written in 
manifestly real form, Eq. (51). However in that form, it is not 
obvious what the generalization i6 to arbitrary m states. The 
generalization is nevertheless simply obtained. Let 
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and 




then the form (51) reduces to that of Eq. (47) for m = 0. For 

r>°'" 

arbitrary m one then need only replace the co £ functions by 

ST'* n k 

the appropriate tOg_ functions, the radial J functions re- 
maining the same. 

Alternatively one can define generalizations of the 


2' 1 
00 C f 


Eq. (49), for arbitrary m. r 

* ex"' o 


with 


3 




t 


v 7 o 


X = "l 

^C h \o) + Cv ’"/° Cy^' 0> ' 

Jt = 3, ■ X * 0 


for m = 0. For X = 0 


The complete function for arbitrary m can then be written 

= +• ^ l ] 

y 

y ^ 

In this case the "radial" functions f^ are the same as in ( 51 ), 

hence real; whereas the angular functions become altered. Note for 

. r\ hr) ± 

/t 0 , however, that the modified spherical harmonics, ^ 


m 


Vt 


are no longer real. 


VI. THE KINETIC ENERGY 

Just as in the case of the angular momentum, the kinetic energy 
can be obtained by a straightforward process of partial differentiation. 


(54a) 


(54b) 


(55) 
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In this case, however, since second partial derivatives are involved, 
the differentiation is a much longer job, and, as we shall see, the 
partial derivatives involving 9 12 no longer cancel out. 

We start then with the kinetic energy in spherical coordinates 

v, 1 = «_ + • i ^ + 1 r i . _l_ £ 1 

r * ^ ^ X* L Sm^ f 

“t _l r i i in •+ i 

r x Ls^ y rt x 

The first two terms are, of course, unaffected by the transformation. 
The angular differentiations then involve the transformation from the 
variables ©i, 9 Z , «=f* 2 to 9, , y , and ©i 2 . 

Consider the coefficient of the ri“ 2 term. After some regrouping, 
we can write 


i_ 

Su I?, 






J_ . ^ 



•WJS; 


+ ^ T ^ j/ 4. — ^ -+ c*<- 3 ** 1 O 

+ i^r r ij* + j- v** lf/»1 


(56) 


(57) 


where for a - 1, 2, 5> 4, refers to 6, ^ , y , ©i 2 . The 
problem thus reduces to finding each of the square brackets separately 
in terms of the Euler angles and ©i 2 . The results are given in Table I. 
The kinetic energy thus becomes 
r 


V, + V, * J 


^ -t JL D_V, + 

r . -*T t n. an 



) i 


12 




( 58 ) 
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where 


F. . 


- — j S*. (Y + ©i» ) 1. 4- Co* (y + ©,» \ Cnt 6 3. 

s,’d ) v I 1 ' >e‘ ' > e 

i 

4 + S.S (t* 4 0,0 Co+9 1_ 

s!b-*V s>" e '»*£ 


(59) 


_ Sw(T«y 4 , 0 ,^) j_ 4 . S*. ( 1 ^ 4 . 8 ,,.) Cot. e ±j 

•Si'hB 'SQ Vp<>6 


and 


- 2 c t# (y + ^ ^ Cot- 0 


A — J_ 


s« 0 


c*t_e c. s (y + e,o 

Sv**e (1 1 


i 1 - i 

J -a<4 


■5 + *■ £ + B. '& 

v >., * 1 " . 




( 60 ) 


ft, 


- Cosjy s~(y + jU - Co ^e So. (i ip + a.*) 


s^e,,. 


So &I*. 


4 Co* 611. ( 1 - J_ \ 

2 e„. Co* e,,. ) 


( 61 ) 


The expressions for F 2 , A 2 and B 2 can he obtained by replacing 0 i2 
by - 0 12 in the above formulae (including the appropriate partial 
derivatives ) : 


F 2 (®> $ > Y * O12 ) =* Fi ( 0 , ©12) 


(02) 


It is clear, since all the coefficients are independent of £ , 


that M coimnutes with the kinetic energy. We have also explicitly 
z 
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verified that Vf + V 2 «J = o. 

Note that the partial derivative involving © 12 and no other 
angles has been placed in the curly brackets with the radial derivatives. 
This is because this term, as the radial derivatives themselves, do 
not affect the orbital angular momentum, and are the only terms which 
act on total S states^’^’'*. 

In fact, in the action of the remaining terms on the angular 
momentum eigenfunction rests the bulk of the reduction of the 
Schrtfdinger equation to its 3 -dimensional "radial" form. With this 
reduction in mind (cf. next section), it is convenient to write F x 
in terms of operators whose effect on the angular momentum eigen- 
functions is particularly simple. One can show 

Fl ■ | Bin M 2 + cos 0 12 (sin 2«f -cos 2w A x ) 

+ sin 0 12 (sin 2t Ai + cos A 2 ) 
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and M 2 i6 the total angular momentum squared operator given in 
Eq. (4l). F 2 is again derivable from Fj. by replacing 0 12 by - 9 i2 . 


VII. THE REDUCED OR RADIAL EQUATIONS, ATOMIC CASE 

The essential properties of the combinations of the operators 
appearing in F x and F 2 , Eq. (62), are the following (cf. Appendix I): 

K 


»r- - ’ 0 - S , SJI-iK) X 


, - . .0 <Jx+V + 

• A/ X. ^ 

D + Wcis-ii] 2 9/ iJr * 1 


A, 4 C„»i^ Aj) ^ 1 ) B/y C< - ) (\ - ^ 

- £.* (• 4 hj, 

*** 


(X-2)- 


{r* 7) - 


(65a) 


(65b) 


W~ 

(s»4 2^ A 1 - C» s iq/ A, ft /jr (l- ^)(| - $• ) 2 


(xV- 

J v t 


, r «. x- 4^ ./ tx+V- 

- V 0- ^ )3r _ a-s,. ,; /\ y /i, +1 ± ^ 

f ['+^(ri-n] 2 V + z 


(S*«<pA ( + 2 if ^ = I*(x-0B, X (»- 


(*-*> 

'/ 


(66a) 


- Si,x (j-2/^ )^ + _ A^ A* +j < 


3 


(J’rl) 


,4 


2 ft 


i 7 


« 2 


(66b) 
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where 

f\ , - C*' O i x f € +i) 

2 ( Mi-0 

and 

- f(f- *+*)(/•*-*)(*+ *-»jj *- 

_____ — 


( 67 ) 


(68) 


Recall that h is the absolute value of k, Eq. (50). We have explicitly 


verified that these relationships are not altered if one replaces the 

<J' r, ' x > ± 

4 1 functions by the J functions of Eq. (, 54 ) .(cf. Appendix 

II) 

This is, of course, necessary for the radial equations to be indepen- 
dent of m. With these relationships, it becomes quite simple to de- 
rive the reduced equations from the original Schrdi dinger equation 


H <P (0 = E V to • 

where the wave function iB expanded in Eq. (51). One obtains 

[L e , a + ce-v)] $r t -/l + i-\r J i- J « 
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L q is the S wave part of the kinetic energy, and only the term con- 
y 12 

12 3 

taining it survives in the description of S states: * * 



L + j. X f x 

t » T % *r t 
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Eq. (70) are the "radial" equations, which it has been our pur- 
pose to derive. They pertain to both types of parity and exchange 
states. Parity is determined by the eveness or oddness of Jt 

If, for example, X is even, and we want to describe a state of 

y-f ( 

even parity, Eqs. (JO) couple the functions and i / for 

jc , o , 1 , U JL Wiis involves l/ x pairs plus one function 

(for X = 0, is zero hence can be taken to be 

zero ) or C + 1 functions . The odd parity equations for the same £ 
correspond to the coupling of the function with x = 1, 3, •••, 

£ - 1. This relates -fy 2 pairs or £ functions to each 

other. Both even and odd parity together therefore involve (2 £ + l) 
functions corresponding the (2 l + l) degeneracy of the vector spherical 
harmonics for a given m. For £ odd, there are -£ functions in- 
volved in the even parity equations and X + 1 functions in the 
odd parity equations. 

For a given parity and l , both singlet and triplet (space 
symmetric and antisymmetric) states are described by the same set of 
equations . The differences in the solutions devolve from the dif- 
ferent boundary conditions which must be applied, Eqs. (52). One of 

v* 

the key virtues of the functions L (ri, r 2 , ©12) is that they 
are either synmetric or antisymmetric; thus they may be confined 

to the region, 6ay ri r 2 . If, for example, t is even so that 

, . 

the exchange character of 1/ is symmetric (which, according to 
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(52), implies that Jy is antisymmetric), then these properties 
may be embodied in the boundary conditions^: 

flL ( ?, ,fi , &,i ) j = o , 

where represents the normal derivative, and 

U r <\ f > . , m ] = o , 

J r,,y 1 

and the solution from there on involves only the region r x y r 2 7 / 0 . 

Such equations have distinct advantages from the point of view of 

12 

numerical solutions 

One can define, however, an asymmetric function in terms of 

which the radial equations can be more simply written. Letting 

* X + *- 

( A • Tii * 6|i) ^ ^ C i Tl ,®i x ) ♦ i N i 

and 

~ * / + x- 

. ft , ) - iy i y, , t-i , e,x ) - e lz ; 

so that, from Eq. (47)# 

fV*,r, ,e,x) •= ± C-0 f/ct, , r, , e»o 


(72) 


( 73 ) 


( 74 ) 


( 75 ) , 


( 76 ) 
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ond obtains by adding (TO®) and ( 70 b) : 

L L - * ’-£ (£ *' • (V * V) [ ) f ' 

- S.,n f'ote,! ( « - 1 A x ) F ( + »i) I F( - (-7— 1 ^ * 
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7*1 


(77) 

e. 
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4.11 . i \ r X & . x (c«jw -c*te„ - j — u/- ^C'- 5 A x)^ 

W V 1 [ ■a*.! 2^ eil 

-v (?-')Jx4 L ; U A ^aL/ { <y * &,„■* C^Hl 

<“ - > *1 i n . _ r* _ / . — ~» 


*' ■* Si H ,tS ..... (E -„]J 


These equations, depending as they do on F { and are more sna- 

2 

logous to the form the P-wave equation of Breit. The question may 
arise in connection with these as well as Breit' s equations, of whether 
they are well-defined, since they involve two functions and IjJ 

and yet there is only one equation (for a given jt ). This question, 
in fact would appear to be particularly relevant as the previous form 
of our equations, ( 70 ) do constitute a coupled 6et for a given X 
To see that both situations are meaningful and in particular that (77) • 

is well-defined, consider a numerical solution of (77) • In that case 

the space of the independent variables is divided into a grid of 

jr 

points, and F^ is the collection of numbers associated with these 



grid points. F g can therefore be considered a vector with as 

many components as there are grid points. The differential equation 
is replaced by a matrix which operates on the vector F* . Now 
everytime an F^ occurs in the equation, it is completely clear 
what has to be done: namely one must let the matrix counterpart of its 
coefficient in the differential equation operate on that component of 
F^ which is its reflection point defined by ( 76 ). This is a com- 
pletely unambiguous prescription which is tantamount to saying that 
the set ( 77 ) is well defined by itself. The reason that (70) is 

composed of two equations for each jc whereas (77) is not is due to 

x 

the fact that the functions are asymmetric and therefore must 

be solved for in the whole Yi, r 2 , 9 12 space. On the other hand 
the f** functions are either symmetric or anitsymmetric, and there- 
fore they are restricted to the ri r 2 , 0 12 (or equivalently to the 
n r 2 , ©i 2 ) space. Since this is only half the independent 
variable space, it is necessary that there be double the number of 
functions to recover the same information. This is again to say 
that (70) and ( 77 ) are completely equivalent. (Nevertheless a re dun- 
dant equation with F^ and F £ interchanged may readily be derived) 
We have stated that (70) has certain advantages from the point 
of view of numerical integration. However, it should also be stated 
that the form ( 77 ) will probably be more advantageous for ordinary 
variational calculations. This is because if one adopts a specific 
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analytic form of , one need only interchange rj. and r 2 in 

the expression to obtain F/*' 

The differences in the f* description from that of 
gives rise to characteristic differences in the formulation of 
boundary conditions for scattering problems (cf. the next section). 

The restriction of these equations to the atomic case (two 
identical particles in a fixed central field) has implicitly been made 
by assuming that the potential is a function of the residual coordinates, 

V = V(rx, r 2 , r 12 ), (78) 

so that V commutes with the angular momentum and therefore appears 
as an additional diagonal term in the radial equations. 

The inter -particle distance ri 2 is related to the independent 
radial coordinates that we have thus far considered, ri, r 2 , © 12 via 
the law of cosines: 

ri 2 2 = ri 2 + r 2 2 - 2rir 2 cos © 12 . 

Alternatively, however, one can consider ri, r 2 , and rj. 2 as the 
independent coordinates and derive radial equations involving them. 

Those coordinates, in fact, have certain advantages since the three l 

singularities in the potential occur at their null, point. As 
such they can describe the wave function in the region of close 
interaction very veil. These variables, therefore, are particularly 
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suited to calculation of low-lying bound states of two-electron 

atoms (where on the whole the electrons are quite close to each other 

and to the nucleus) and such successful calculations have been done 

ever since the early work of Hylleraas . ^ 

When one considers the equation in the form we have previously 

given them, involving ©12, one is naturally led to expand the "radial" 

wave function in terms of Legendre polynomials of cos . The 

expansion is then truncated after some P n (cos ©i 2 ) and convergence is 

sought as a function of n. In these classes of two electron problems, 

this constitutes the idea of configuration interaction in its most 

13 14 15 

general form. Recently this idea has come under some criticism * * 

principally because such a relative partial wave expansion necessarily 

converges slowly where the electron-electron interaction is large 

(ri 2 small). The argument is doubtless Justified for the above- 

mentioned low-lying bound states. However the argument can easily 

get distorted and exaggerated, for instance when applied to the 

14 

low-energy scattering of electrons from hydrogen. The point there 

is that the long-range correlation coming from the induced potential 

in the atom is at least as important as the short range correlations ^ 

and yet is only poorly approximated by the conventional Hylleraas 

12 

type of expansion. This situation has been discussed elsehwere . 

These reservations notwithstanding, however, it is nevertheless 

true that the most accurate three body calculations have been made 

17 

using the ri, r 2 , r 12 coordinates, or linear combinations of them , 
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lo '* 17 

on- the low-lying states of helium and its isoelectronic ions . 

We therefore give below the radial equations in terms of ri, rx 2 * 

The equations are in their asymmetric form corresponding to Eq. (77)* 

since it is assumed that they will be utilized in connection with 

variational calculations with analytic expansions of the radial wave 

functions . 
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Here 


% - ♦ [^i )* * n a ;J * 


(80) 


The quantity whose square root ^ is can easily be shown to be 


pc 

positive definite. In the equation (68) the iB understood 

to be a function of r lf r 2 , r i2 : 

K - F /(v ^ 

In addition is the kinetic energy counterpart of the 

S-wave in terms of r x , r 2 , r i2 : 


L = l + i 
^ r. 


* £ T tl 

t, n : . 


.2 

'fa. 


4 - A . A. 

* \ 




+ £ 

^ r, ^ 


(81) 


The equations ( 79 ) can readily be put in the form of coupled 
equations for a given / . In that form they would be closest to 

the form originally given by Holmberg^, ( although as we have stated 
his equations apply to three equal mass particles). One salient dif- 
ference between the two sets of equations, however, is that the 
present ones are manifestly real, whereas one term in Holniberg's 
equations is imaginary . It is clear that the equations as well 
as the solutions must be reducible to completely real form for any ♦ 

given angular momentum state. The accomplishment of this in the 
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present case comes from the explicit construction of real vector 
spherical harmonics, Eq. ( 49 ). 

VIII BOUNDARY CONDITIONS FOR SCATTERING 

In this section we derive the asymptotic forms of the 

radial functions corresponding to the scattering of an electron 
from a one electron atom in its ground state. The Coulomb modifica- 
tions when the target system is an ion instead of an atom can 
readily be made and will have no effect on the angular Integrations 
with which we are here concerned. 

As we have seen in the foregoing sections culminating in the 
last section, the selection of a symmetric choice of Euler angles 
(Hdriberg's angles) has allowed for a completely general derivation 
of the radial equations. From the point of view of a scattering 
problem, however, a symmetric choice of angles is not the most ad- 
vantageous since here we are concerned with an intrinsically asym- 
metric situation. Thus if we consider that region of configuration 
space where ri is large and r 2 small, corresponding to electron 1 
being scattered from the atom to which electron 2 is bound, the 
wave function in this region alone will not be symmetric. However 
in terms of the Hylleraas-Breit angles, the spherical angles of one 
of the particles being defined as two of the Euler angles, the wave 
function in this asymmetric region is easier to describe. Neverthe- 
less this is a complication of detail only, since all the angular 
integrations may readily be performed as we shall now show. 
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We start with the statement that the complete wave function 
must have the asymptotic form: 

Lvtn V i / £o ( £ , ** ) ^ y t0 (St, ) 

n £ 

where (Ri B (r 2 )/r 2 ) Y 0o (&») is the ground state of the one -electron 
atom (hydrogen). On the other hand, from Eq. (51) » 


(82) 


* X+ X+- 

Li- = s ^( K/ » + *< ~ i] ?) 2. [•*, (*„)&(»,£ M 




where 


r* * 

X - 


(83) 


+ «(* (e,x) & >*)J 


m 5 %* 

e<^e at a^ 


( t?ix) - X> 


(84) 


It should be noted that (82) refers to the state of parity (-1) 
as long as we are considering elastic scattering from the ground (is) 
Btate. ThiB then defines the evenness or oddness of the values over 
which x goes in the summation in Eq. (83). 

The quadrature in (84) can readily be performed by recalling 
from Section II that 2^ is the angle between £ and r L 
whose spherical angles in the primed coordinate systems are given 
in Eq. (20). One can then use these spherical angles to expand 


(cos ft) via the addition theorems for spherical harmonics. 
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In its real form^ this gives in the present case 

l m 

v ( W - P(T i’J'j.) fpve) + l i .(-0 ( L'7?±' <^ a ) P?7e) ( nuif- J-6, a )) 

■">=« (O m)l X ' 


( 85 ) 


In ( 85 ) we have written both the Legendre and associated 
Legendre polynomials as functions of the angle but what we mean in 
all cases is that the angle is to be substituted into the trans - 
cendental form of the function. For example Pi (p ) = cos p> and 
not Pi ( p ) = p . The sign of the P c m is that of Magnus and 
Oberhettinger^ (which differs by (-l) m from that of Morse and 
Feshbach^). To complete the quadrature in (84) we note that 
o,* 


oL, ie) - c-0 l( t-*u V p £ (e) 
L ft +x)i J 


(86) 


Substitution into (84) now yields the desired result: 

4(V) * [ r i^| y [ + 0 - co» ^ ^ 

S- <i»,a J 

*■ L (l+X)\ ( J 


(87a) 


( 87 b) 


The radial function themselves thus approach 

TCrt 


f ( ( T, , t, , e ri ) = S-Ut! ) R U (V, 

n ~ ' ,v) 


( 88 ) 



in which form we see that ri, r 2 dependence of all the 
functions is independent of / , so that none of them vanishes in 

the asymptotic region. Since in all cases the 612 dependence is 
trivial, it may be worthwhile to define new functions whose asymptotic 
behaviour is strictly the rj., r 2 dependence in (88). 

For bound state problems, it is clear that all the radial func- 
tions must vanish in all asymptotic regions. 

IX OTHER APPLICATIONS 

In addition to two electron atomic or ionic systems the present 
equations apply to double mu or pi mesic atoms, although as the mass 
of the identical particles gets heavier, the correction for the center 
of mass becomes more important. Also for the spinless bosons (pi 
mesons) only the space symmetric solutions will presumably be relevant. 

The equations can also be applied to two different particles 
of the same mass (positron -hydrogen scattering, for example). In 
this case, the potential, V, will no longer be symmetric hence the 
solutions will not be symmetric which implies that boundary condi- 
tions like (72) must be changed to matching conditions of the asym- 

19 

metric solutions along the line ri = r 2 . This has the effect of 
giving one solution where formerly there were two, in accord with the 
distinguishability of the particles. 

The major extension of thiB approach is to two-electron diatomic 



molecules. In this case, the extension from one to two elec- 
trons is non-trivial. However, the analysis has been completed 

21 

and will be published elsewhere. 
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Appendix I 

In this appendix, we prove the Eqs. ( 65 ) and (66). For in = 0 
^ terms give zero. Therefore, we can write 

A, — > ^ _ Cot 0 "sL — Cot 

*e 2 a© 

_> 1 Cot 0 ?£ - L _ 2 Cot - . 5. 

bipH ** '*1' 

and 

x + X 

3/ = N, Co>x<f 

S^ICiy Si-B Fy 

where 

n< x s /zg+t (775 7" . i_ j_ 

J *▼* (* -/)/ *! 

and 
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F x = F , x + * +1 , I «-* , So!© ) 
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(15) * 


( 16 ) 
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Now 


A, 3, = n„ CosXtf/^ ix{ jl-\) Sv^e Cos8 f> - 


x s«i6 f* 


+ 1 

X J r * k e J c , 

x <mGC os© by r x ^ ( +■ Su.© Ay ^•f-2j 

(17) 

where we have used the well known relations for the derivatives of 
hypergeometric functions: 


if, „ A 1 

de 

St*0 Pjr -t* | 

( 18 ) 

d>, . a‘ 

do 1 

c «8 F x+( + A x Sw© F x ^ 

(19) 

A relation between F x , 

F Xl _ ( and F x +z can be obtained from 



the differential equation satisfied by the hypergeometric function 


C ~ e fxti = F x - 


a { *1" 

* ^4-2 
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Using (110) in (17), ve find 

Xf 
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O'. 


Similarly, . 


- ” N* Sv- x«f 




X 

x S-ce 


-t- 



- 

X+l 


Xti - 
.Svw 




( 112 ) 


Multiply (112) by SW2f and (111) by cos 2t and subtract to get 


( t \^ X+ e t *{** l) + 

( Sv* /\ a - c-» 5^ A, ) 2 f = - Co>(*-2)Y q - Nt* A, A jr+I 4, (113) 

lK W 
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Let x — y r + I 


in (110), then 


Co S 0 


- F- 


- S«s© F. 




Substituting in the above for F >+ , 
we get after some rearrangement 
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e A 1 

a(* +/)■*- 


and F 


X*3 


by using (110) 


Sv« © ft 
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Letting % multiplying by 2( X — I ) 2 sin 2 © and rearranging we 

have for y / ‘}_ 




G 
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^0,1 

(115) 

Also we find 

directly from (111) and (112) 



G 

= i e hi F a 

)f C O 

(116) 

G 

= - (j A*) Sue F, 

*: - I 

(117) 


Finally then with the substitution of the above in (113) we. obtain 


for K&Z, 
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x O'-*) 

(S» 2^ K i *(*-') ft ty 


Where Q/*' has already been defined in Eq. (68). 

The special cases x - o, | can be determined from (113), 

(116), (117)* With proper normalization 

(S«* iV'A.- c -»^A 1 )X 4 = -/I 

and 

(S^tO'Aa - Co%iy = (j-iA* J^ + - Af Al ^ 

We can combine (Il 8 , 119,1 £0) to get Eq. (65a). 

Similarly, we can prove (65a), (66a), (66b). 
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Appendix II 

For m ^ 0 , we get, using Eqs. ( 44 ) and ( 64 ) 
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For illustration purposes, let . Multiply (112) by 

sin 2f and (ill) by cos and subtract to get 
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Using fris result, (which reduces to (115) for m = 0) in (115), we 
obtain for k £ 0,1 
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It is easily verified that 
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Using (1110) and (1111) in (Il8), ve find 
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Similarly, ve find 
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Adding (1112) and (1113) and using the definition (54a) ve derive 

finally for * ± 0, 1 
\ 



( 1114 ) 






c\ " **■ 

feiif A,- C*^A, W, »- a x^r-oB # 3^ - AjjC, 
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’Or*a 


which is identical in form to (118) for m = 0. Other relations can 
be proved in the same way, thus the radial equations are in fact 
independent of m. 
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Table 1 : Coefficients of the Angular Derivatives in the Kinetic 

Energy a Cf. Eq. (49) 


Coefficient 

Derivative 

Coefficient Derivative 

sin 2 (Y + — ) 

a 2 

cot ©12 

a 

2 

sin 2 ©i 2 
cos 2 (y + 

sin 2 ©i 2 sin ^ 

a©* 
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a©i 2 
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Ai 
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sin 2 ©i 2 

aeaY 
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d© . 

-2 cos © cos 2 (y+^ 2 ) 
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O 

a 2 

sin 2 © sin 2 ©12 
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V 

a$a©i 2 

Bi 

a_ 

-1 

a! 


a* 


3¥3©i2 


a Ai and Bi are given in Eqs. (60) and (6l) 


Figure Caption 
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Figure 1. Perspective drawing of (Holmberg's) Euler angles 
and the unit vectors of the problem. 
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